Anomalous Absorption Line in the Magneto-Optical Response of Graphene 
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The intensity as well as position in energy of the absorption lines in the infrared conductivity of 
graphene, both exhibit features that are directly related to the Dirac nature of its quasiparticles. We 
show that the evolution of the pattern of absorption lines as the chemical potential is varied encodes 
the information about the presence of the anomalous lowest Landau level. The first absorption line 
related to this level always appears with full intensity or is entirely missing, while all other lines 
disappear in two steps. We demonstrate that if a gap develops, the main absorption line splits into 
two provided that the chemical potential is greater than or equal to the gap. 

PACS numbers: 78.20.Ls, 71.70.Di, 81.05.Uw 
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While graphene has only recently been isolated 1 -, 
several anomalous dc properties already appear 
well established 2 -^. This includes an unconven- 
tional integer quantum Hall effect with filling fac- 
tor v — ±2, ±6, ±10, . . i£ i 3 i 5 i 6 i 7 which changes to 
v = 0,±1,±2,±4, ... at a magnetic field B above 20 T 
Ref4, which may indicate an opening of a gap. Also 
observed is a Berry phase of ir—£& which increases to 
2ir in bilayers 8 . These spectacular dc properties can be 
traced to the nature of the quasiparticles in graphene 
which are massless Dirac fermions governed by 2 + 1 
dimensional quantum electrodynamic o 9 ' 10 . Optical 
conductivity measurements play a central role in solid 
state physics and have provided additional information 
on the dynamics of charge carriers and its change 
under external perturbations, not directly available in 
dc measurements. Here we present a theory of the 
diagonal optical conductivity of Dirac quasiparticles in 
an external magnetic field and highlight several unusual 
features, not present in a conventional two dimensional 
electron gas, which bear directly on their unconventional 
dynamics. We also investigate an emerging new physics 
of graphene related to the sublattice symmetry breaking 
reported in RefA 

In the presence of a magnetic field, Landau levels (LLs) 
form in the electronic density of states and transitions 
between these levels give rise to absorption lines in the 
optical conductivity. In this letter we find that the se- 
lection rules, which set the frequencies of the lines and 
their relative spectral weight, provide a distinct signature 
of the presence of the Dirac quasiparticles. Not only do 
the positions and the intensities of the absorption lines 
scale with \f~B~ rather than B, but the line corresponding 
to the transition from the lowest, n — to the n = 1 
(Landau level) LL is anomalous due to the special Dirac 
character of the n = level. This level and no others can 
participate in both inter and intraband transitions. As 
a consequence, the first absorption line always appears 
with full intensity or is entirely missing, while all other 
lines can also be seen to have in addition half intensity 
depending on the position of chemical potential \i relative 
to the LL energies. The chemical potential fi is tunable 



in experiments by application of a gate bias voltage to a 
field effect devicei*2i2il. 

Expressions for the absorptive part of the diagonal as 
well as Hall conductivity on which our calculations are 
based can be found ir* L n i 12 . In the limit of zero impurity 
scattering rate, r — > 0, the diagonal conductivity as a 
function of the photon energy Q is 
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where the thermal factors are A(T) = n F (—M n +i) — 
n F {M n ) + n F (-M n )-n F (M n+1 ) and B{T) = n F (M n )- 
n F (M n +i) + n F (—M n +x) — n F (—M n ) with the Fermi dis- 
tribution n F (uj) = [cxp[(cj — fi)/T] + (we set the 
Boltzmann constant fc^ = 1). To include damping, the 
n'th delta function in Eq. |T]) is replaced by a Lorentzian 
of width r n + r„+i with r„ the width of n'th level. Also 
in Eq. |T]) a term fl — ► —17 is to be subtracted. In the 
presence of an excitonic gap A 13 i 14 the LL ene rgies are 
E n = ±M n , with M n = ^A 2 + 2n\eB\hv F /c, and v F 
is the Fermi velocity in graphene, h is Planck's constant 
h over 2-7T, c is the velocity of light and e is the charge 
of the electron. The field B is applied perpendicular to 
the plane and for B = 1 T the energy scale Mi w 420 K 
for v F « 10 6 m/s and A = 0. As for the conventional 
case, the selection rules allow transitions only between 
adjacent LLs (see Reffii for details). 

Our numerical results for Re <j xx (O) are given in 
Fig. D](a), where it is plotted in units e 2 /h as a function 
of £1 for A = 0. The same constant electronic scatter- 
ing rate r = 15K 4 was used for all LLs, an assumption 
that could easily be relaxed. Except for the long dashed 
(red) curve for which the magnetic field B = 10~ 4 T, and 
is shown only for comparison, B = IT for other curves. 
They differ only in value of chemical potential [i. The 
dash-dotted (black) curve has M < fx < M\ ; solid (blue) 
curve has Mi < fi < M2, and the short dashed (green) 
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FIG. 1: The absorption peaks and corresponding transitions. 

(a) Real part of the conductivity, Re a xx (Q.) in units of e 2 /h vs 
frequency Q in cm -1 for temperature T = 10K and scattering 
rate F = 15K. Long dashed (red), the chemical potential 
H = 50K and the magnetic field B = 10~ 4 T, dash-dotted 
(black) jj, — 50K, solid (blue) /i = 510K, short dashed (green) 
fi = 660K all three for B = IT. (b) A schematic of the allowed 
transitions between LLs n = 0, 1, ... 4 shown as heavy (red) 
dots for both positive and negative Dirac cones. The pair of 
cones at points K and K' in the Brillouin zone (see Fig. 3 

(b) ) are combined. Three values of chemical potential are 
shown as the (violet) heavy solid line. At the left, fi falls 
between Mo and Ml; middle, between Mi and M2; and on 
the right, between M2 and M3. The first line on the left 
which corresponds to the first peak in dash-dotted line in 
upper panel is unusual in that it arises from both interband 
and intraband transitions and ceases to occur for all other 
values of chemical potential, while all other interband lines 
first drop to half their intensity before disappearing entirely 
as fi crosses to higher energy levels. 



sitions from negative to positive Dirac cones. As fj, is in- 
creased to cross the n = 1 LL [solid (blue) curve] , the first 
peak in the dash-dotted (black) curve disappears and the 
intensity of the next higher peak drops to half its value, 
while all other peaks remain the same. In addition, a 
new intraband peak appears at Mi — M\ and its inten- 
sity picks up the lost optical spectral weight from the 
missing and/or reduced interband transitions. Compar- 
ing with the short-dashed (green) curve, the first two in- 
terband lines have disappeared, the third has halved and 
the intraband transition is at M3 — Mi. The first line of 
the interband series is different in that it never appears 
with half intensity, while the others do. It is the only 
one which can be considered as involving both interband 
and intraband transitions as can be seen in Fig. [T] (b). 
The energy level schematic shown in this figure summa- 
rizes our results. It shows how the first interband line on 
the far left appears alone, while all others come in pairs. 
This transition exists only for values of the chemical po- 
tential less than Mi, while the others first drop to half 
intensity as another LL is crossed before disappearing 
at the next crossing. To verify these predictions optical 
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FIG. 2: Real part of the longitudinal conductivity, Rea xx (Q) 
in units of e 2 /h multiplied by F / \J heBvp / c as a function of 
normalized photon energy Q / ^2heBvp/c. The temperature 
is T = 10K, the scattering rate F = 15K and the chemical 
potential fi = 680K. The long dashed (red) curve is for B = 
0.75T, dash-dotted (black) for B = IT, solid (blue) for B = 
1.5T and short-dashed (green) for B — 3T. 



curve has M2 < \x < M3. The curves remain essentially 
unchanged for any values of fi between Mn and Mjv + i 
and evolve rapidly from one pattern to the next as a new 
LL is crossed, for a change in chemical potential of order 
T or r whichever is greater. 

Referring to Fig.QJa), m t ne dash-dotted (black) curve 
the lowest energy peak is at M\ and the others at 
Mi + M2, M2 + Ms, ... all correspond to interband tran- 



experiments on graphene based field effect device would 
allow one to change the chemical potential as in Fig. [1] 
through a change in gate voltage. Such optical experi- 
ments have been reported by Li el alJ^ in organic metals. 
It is however also possible to get the same information 
for a fixed chemical potential if Hea xx (Tl) is measured 
at a few well chosen values of the external magnetic field 
B. This is illustrated in Fig. [5] To see the desired effect 
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it was convenient to normalize to y // 2heBvp/c and to 
multiply the vertical scale by the dimensionless quan- 
tity T/^heBv 2 F /c. This arrangement guarantees that 
the various interband lines will all fall at the same nor- 
malized value of Q/ y/2heBvp/c and the intensity of the 
various lines follow the same pattern as in Fig. [TJ 

To understand the intensity that is to be assigned to 
the various lines and their changes with (i it is sufficient 
to take the limit T — > in Eq. JTJ) . For chemical potential 
/z € [M N ,M N+1 ] we find (ft > 0) 
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where the first and second terms give the interband and 
intraband transitions, respectively, at the energies indi- 
cated by the 8- functions. The denominators come from 
the I/O factor in Eq. ([J). In units of (e 2 /h)M 1 (n/2) 



which goes like V B, the optical spectral weight of various 
lines is given by the quantities multiplying ^-functions in 
Eq. |(3J|. For N > all interband lines below n = N are 
missing and the intensity of the line at n = N is half of 
2/(yn + 1 + y/n). The total optical spectral weight lost 
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The right hand side is exactly the intensity of the corre- 
sponding intraband line given by the second term in {2} . 
We verified that the integrated optical spectral weight 
under Re a xx (f2) up to energy Q m ax is linear in O maa; and 
approximately equal to (e 2 / 'h)(ir /2)Q max provided £l max 
is taken large enough for many LLs to be involved. This 
linear in £l m ax behavior which holds as well for B = is 
another signature of the existence of the Dirac cones. In 
ordinary metals for B = the optical sum would rapidly 
saturate to a constant value as £l m ax is increased beyond 
the scattering rate T. The above rules for the intensity 
as well as the peculiar behavior of the first inter- (intra- 
) band line as fj, is increased from zero are specific to 
graphene and as yet have not been verified experimen- 
tally, although some data exist as we now describe. The 
early magneto-reflectance data^ in fields of order of 1 
to 10 T in graphite assigned LL lines to the ii-point in 
the Brillouin zone. Some of these have central frequen- 
cies which vary as the square rooti£ of B as expected in 
graphene. More recent data by Li et al<^ who use high 
fields up to 20 T, show conventional linear in B behav- 
ior. However, very recently, infrared transmission data^ 
on ultrathin epitaxial graphite samples^ in fields up to 
4T have become available and provide Kea xx (fl). A \[B 





FIG. 3: The absorption peaks and the corresponding transi- 
tions when the lowest LL is split, (a) Real part of the longi- 
tudinal conductivity, Rea xx (Q) in units of e 2 /h vs frequency 
n in cm" 1 for temperature T = 10K, L = 15K, B = IT 
and chemical potential jj, — 150K for 4 values of the exci- 
tonic gap A. Long dashed (red) A = OK, dash-dotted (black) 
A = 100K, solid (blue) A = 150K, short dashed (green) 
A = 200K. (b) Possible optical transitions between LLs with 
an excitonic gap A. Two pairs of cones at K and K' points in 
the Brillouin zone are drawn. Only LLs with n = with en- 
ergies —A (at K point) and A (at K' point), and n — 1 with 
energies ±A/i , respectively, are shown. Three values of chem- 
ical potential are considered \fi\ < A (green), /i = A (blue) 
and > A (black). For the case /i = A the state at energy 
Eo = A is both occupied and unoccupied with probability 
1/2. 



scaling for both the energy and the intensity of the first 
line was taken as evidence for Dirac quasiparticles. In 
view of the known sensitivity of some properties to cou- 
pling between layers it is not obvious that the various 
carbon sheets in the samples are sufficiently decoupled to 
be treated as isolated graphene sheets^. Further, to un- 
derstand the results it is necessary to assume that there is 
a variation in charging, i.e. each sheet can have different 
chemical potential. We have analyzed the intensities of 
the one intraband and of the three interband transitions 
reported m}~. The intensity of the first interband line (B) 
in Fig. 1 of Ref. 18 , which can occur only for < Mi, is 
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the largest and that of the intraband line (A) is smaller 
by a factor of 0.36. This can only be understood if there 
are at least three times more sheets for which the chem- 
ical potential falls below Mi than above. Further, the 
ratio of the intensity of the second (C) and the third (D) 
interband lines to the first (B) are 0.34 and 0.19, respec- 
tively. These values are somewhat smaller than 0.41 and 
0.32 expected from our theory, coming from the sheets 
with < Mi. Adding in the other sheets with > Mi 
only increase the above discrepancy, because they add to 
the second and third interband lines, but not to the first. 
While the above provides some indirect support of our 
theory, further experiments are needed to verify the main 
features described here. A case for further studies is a 
bilayer— where it is known& that the two lowest LLs both 
lie at zero energy and this should change the pattern of 
behavior of the absorption lines described here. 

The ac conductivity may also be used to get informa- 
tion on the opening of a possible gap related to the sub- 
lattice symmetry breaking in graphene^ i 13 ' 14 i 22 ' 23 i 24 . In 
Fig- El (a) we show results for Hea xx ($}) plotted using the 
generalization of Eq. |T]) to finite T. Four values of the ex- 
citonic gap A are considered. For A < the single peak 
at 294 cm -1 has split into two of nearly equal intensity, 
for \fi\ = A the lower peak has approximately 1/3 the in- 
tensity of the upper peak and for A > a single peak is 
seen at the higher energy. Fig. [3] (b) gives a schematic of 
the levels involved and shows the transitions. For conve- 
nience we have fixed the gap in the figure and considered 
three values of chemical potential. For < A the inter- 
band transition — Mq to Mi (LHS cone) and —Mi to Mq 
(RHS cone) are both possible and have the same energy 
Mi + A, so only one line is seen shifted to higher energy. 
For |/u| = A one needs to account for the fact that the 
state at n = (RHS cone) is occupied (unoccupied) with 
probability 1/2. The possible interband transitions are 
-M to Mi (LHS cone) with energy Mi + A and -Mi to 



Mo with energy Mi + A, and an intraband transition of 
M to Mi with energy Mi — A (RHS cone) each weighted 
by 1/2 rather than by 1. This gives a 3 to 1 ratio favoring 
the higher energy line. Finally, for > A, the —Mo to 
Mi interband transition (LHS) has energy Mi + A and 
the intraband transition of Mo to Mi (RHS) has energy 
Mi — A. The two lines have nearly the same intensity. It 
is clear that optics provides a useful probe for observing 
various kinds of gaps and here the excitonic gap is taken 
as an example. 

In summary, we have found an anomaly in the first 
absorption peak in the optical conductivity of graphcne. 
When the chemical potential, tunable by the gate volt- 
age, sweeps through the various LL energies, this peak 
never halves its intensity before disappearing while the 
other lines do. This pattern of behavior is specific to the 
Dirac nature of quasiparticles in graphene and should be 
easily observable in the experiments. It originates from 
the special structure of LLs in graphene, where each pos- 
itive energy level has its negative counterpart, except for 
the zero energy lowest LL. For the Shubnikov de Haas 
oscillations^ this structure of LLs resulted in a phase 
shift of 7r used to identify these quasiparticles in the 
dc measurements^- 3 . Furthermore, should an excitonic 
gap open in high magnetic field, its presence will lead 
to a specific signature in an ac optical experiment. For 
A < \fi\ < Mi, the n = to n = 1 transition will split 
into two lines of nearly equal intensity, while for = A 
the intensity ratio of higher to lower energy line will be 
about 3. For < A the line will simply shift upward 
without splitting. 
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